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ABSTRACT 

In recent studies it has become increasingly clear that presence of infinitely many 
instability bands of the parametric resonance plays crucial roles in the phenomenon 
of particle production under periodic classical field oscillation. We extend previous 
works to a general class of models including both the Yukawa and the quartic type 
of couplings of the classical field to quantum bose fields. Decay rate from the n— th 
band is derived in the small amplitude limit using the functional Schrodinger picture. 
It is then shown that this analytic result of the decay rate can also be derived as the 
zero momentum limit of a physical process, n particles that comprise the classical 
homogeneous field decaying simultaneously into 2 bose particles. The latter approach 
uses ordinary perturbation theory, hence the former result is a novel resummation 
of many perturbative amplitudes, which usually becomes complicated for a large n 
order. 



*To appear in the Proceedings of the Symposium on "Frontiers in Quantum Field Theory" , 
(World Scientific, Singapore, 1996). 
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1 Introduction 

Particle production under periodic perturbation is a physical process that 
governs many problems in different areas of physics. In cosmology in particular, it 
is a dominant process of entropy generation after inflation, since inflation dilutes 
away essentially everything in the observable part of our universe and leaves behind 
the field (inflaton) oscillation As another example one may cite presence of 
flat field directions in supersymmetirc models. It may either give good effects such 
as the Affleck-Dine scenario of baryogenesis 0, or bad effects such as the Polonyi 
or the modular problem [B| that may potentially destroy the successful result of 
nucleosynthesis. In all these problems effects of the parametric resonance are crucial 
if the initial oscillation amplitude is large enough. Thus a deep understanding of the 
basic process is inevitable in any of these applications. 

In our previous study 1, || we proposed to formulate the problem of particle 
production by taking a short time average of the quantum density matrix over a few 
oscillation periods. In our view this time average is a substitute for a more physical 
means of the coarse graining when created particles interact among themselves or 
with some other particles. A large fluctuation of the quantum system coupled to the 
classical periodic oscillation is essential to this viewpoint, and the quantum system 
exhibits a nearly classical behavior even before taking the time average. 

In this work we bridge a gap by offering in the language of ordinary perturbation 

of interacting quantum field theory a simple understanding of the decay formula 

previously derived. We do this by generalizing the coupling of oscillating field to 

quantum bose fields beyond that previously analyzed. The model Lagrangian for 

the oscillator coupling is thus taken as 

I i 

with £ the classical oscillating field that takes a simple sinusoidal form, 

= £ sin(m^) . (2) 

is the mass of the £ field, and <p is a generic bose field treated here as a quantum 
field. In the rest of discussion it is not important to assume a particular relation 
between the two dimensionless couplings, the quartic coupling g± and the Yukawa 
coupling gy, but in order to organize results systematically we assume these to be 
of the same order of magnitude; g±/ gy = 0[1\. 
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The quantum state coupled to the external £ oscillation evolves according to the 
functional Schrodinger equation, which is much simplified to an independent set of 
equations in Fourier modes due to the translational invariance of the background 
oscillation. Thus the state vector in each k mode obeys the Schrodinger equation 
with variable frequency, 

u&t) = k 2 + gie(t) - g Y m^(t) • (3) 

Precise relation between the quantum wave function ipk{c[k,t) and the classcial os- 
cillator equation has been established It takes the most definitive form if one 
starts with an initial quantum state of the ground state of some reference frequency 
ipk(Qk 1 1) = (uj/tc) 1 / 4 e~ exp[— | uql] , for t ~ 0. At any finite time t 

I 7T id 

Mlk,t) = — exp[--j — — $ + - — ln|u fc (t)| 2 -^], (4) 

J\Uk{t)\ z \ u k{ l )\ 

where the time dependent complex function Uk(t) obeys the classical oscillator equa- 
tion with the definite initial condition, 

-^± + l j 2 (t)u k = 0, u k (0) = (u/n)- 1 / 2 , u k (0) = iuu k (0). (5) 

The meaning of the classical oscillator amplitude Ufc(i) is thus unambiguous: its 
modulus \uk\ governs the Gaussian width of the wave function and at the same 
time its logarithmic derivative gives the fluctuation of the state in terms of the 
phase factor. This can be seen most clearly in the Fock space base of harmonic 
oscillator of the reference frequency: the diagonal density matrix element derived 
from this wave function has a smallest fluctuation and is related to the average 
particle number (N u ) according to P2n.2n e~ n ^ N ^ / y ivn (Nj) @, while the 
off-diagonal density matrix elements contain wildly varying signs of ±, giving zero 
after taking the time average over a few oscillation periods. This phenomenon occurs 
in infinitely many band regions of the parameter space (k , £o) within which a generic 
classical Uk exponentially grows; Uk(t) — > e Am?i//2 x (periodic function) with A > 0. 
The exponential growth implies rapid excitation of high harmonic oscillator levels. 
With the short time average it can be interpreted that particle production takes 
place with (N u ) oc e Am «* . 

After coarse graining of the time average one can discuss the decay law of the 
initially prepared ground state: it follows the exponential form, 

poo-e- r ^, r = fr n , r n = ^ £ a s . (6) 

n — 1 k en-band 
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r is the total decay rate per unit volume and per unit time. Computation and 
interpretation of the decay rate T n of the n— th band is our main task in the rest of 
discussion. 

2 Small Amplitude Analysis Revisited 

It is customary to recast the classcial equation into a dimensionless form, 
d 2 u 

— + [h-26 4 cos(Az) - 26 Y sm(2z) ]u = , (7) 

z = m(t/2 , , ! = 4 ^! + 2»l#, ft = Mf, » y = 2^, (8) 

ra| ra| m| 

with m the mass of the bose quantum field (p. Following the general theorem for 
solutions of differential equation with periodic coefficients, one expands the solution 
in the form M, 

oo 

u(z)= J2 c k e^ x+in+2ik > , (9) 

k=— oo 

with n = 1 , 2 , 3 , • • • . The condition of the existence of non-trivial solution is 
then equivalent to the condition of non-vanishing matrix determinant of infinite 
dimensions, with the diagonal entry of 7& = h + (A + in + 2ik) 2 , and the next off- 
diagonal entry of ±i6y and the next-to-next off-diagonal entry of — #4 and all others 
= 0. 

In the small amplitude limit of |£o| 1 the n— th instability band starts at h — 
n 2 . The boundary curve h = hS n '(^o/m^) dividing the stable and the unstable bands 
corresponds to the eigen-modes of the form, cos(n2;) and sin(nz), in the £0 — > limit. 
With the assumption that the coefficients, Co and c_ n , are dominant and the rest of 
Cfc, A, and h — n 2 are all small, one may simplify the structure of this matrix such that 
7^ = 7_„_fc = — 4:k(n+k) (k 7^ , —n) and 70 = h—n 2 +2in\ , 7_ n = h—n 2 —2in\ . 
Dividing the matrix into a few parts, one first solves the top (k > 1) and the down 
(k < —n — 1) infinite dimensional parts in favor of c\ , C2 and c_„_i , c_„_2, 

c 2 = D 2 ^(e l c + e 2 c_i) + D 22 l e 2 c , c x = D^{e x CQ + e 2 c_i) + D~ 2 l e 2 co , 

c -n-l = -^11 ( e l c — n "I" e 2 c — n+l) "I" -^12 e 2 c -n ) 

c_ n _ 2 = D 2 l*{e\c- n + e2C_ n+ i) + D^^C- n , (10) 

with ei = — i^y and e 2 = ^4 and the matrix inverse D~ x defined as a limit of big 
matrix in the left-upper and the right-down (identical) corners. One next solves the 
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central block (—1 > k > —n + 1) for c_i , c_ 2 , c_ n+2 , c_ n+ i in terms of c , c_ n . Here 
one needs to invert, ignoring sub leading terms, 



7-1 -ei -e 2 

-e\ 7_2 -ei -£2 

e 2 











o ^ / 








e! 7-n+2 - ei 
e 2 -£i 7-n+i / 



C-l 
C-2 



c 




C-n+2 





e 2C_ n 

\ eic_ n J 



(11) 



One finally solves the determinental condition for the equation cq , c_ n , which 
reads to the leading order as 



det 



h-n 2 + 2mA - A 



C 







.4 



-C* h-n 2 - 2mA - A 

NW + ^n) , C7 = e 1 (e 1 E^ n _ 1 + e 2 E^ n _ 2 ) + e 2 {e x E^ + t 2 E 2 ) n _ 2 ) 



(12) 



(13) 



where E is the matrix in the left hand side of eq.|TT[ 

We regard e\ = 0[e) and e 2 = 0[e 2 } since g 4 / gy = 0[1], and work out matirix 
elements to leading orders of e. The net result may be summarized as the formula 
for the growth rate A, 



A 



n ~ 2n\ 



A 2 n -(h-n 2 



3. )2 A 

2 1 \ J ' 
1 ^ 



IC 



C n = det 



^ -%Q Y 4 

- 7_i - i9 Y 4 







V 



2 2(n-l) [( n - l)!] 2 

N 
••• 




-7-n+2 -iOy &A 




(14) 



(15) 



-7_ n+ l -IU Y J 

with 7_fc = 4&(n — A;). This result generalizes the previous one in ref [f|] in which 
6*4 = was assumed, hence C n = (—i9 Y ) n . The decay rate of the n— th band is 
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computed by summing modes within the band in the narrow width approximation, 

, 1 



m 



256vr \ 



4m 2 



A 2 

1 — *n • 



n 2 m 2 n 



(16) 



The case of n = 1 is an exception to this general formula that must be treated 
separately, resulting in 



Ai(Mo) 



\ U (h i + ^v r - gW & 



1 



4m 2 



;i7) 



3 Physical Interpretation in terms of Familiar Perturbation Theory 

We shall now offer interpretation of the decay rate T n of n— th band in terms 
of the ordinary perturbation theory. Ordinary perturbation uses the particle picture 
and for that purpose it is useful to recast the rate formula in terms of reaction 
rates of indivisual particles by dividing some powers of the particle number density 

= |m^Q. Let us first note the dependence of the rate T n on the oscillation 
amplitude, oc £o™, and the momentum of final ip particles in the zero momentum 
limit of E% — ► rri£ for the process, n £ — > <f<f] p (/3 = ^y 
appears in the rate formula, eq.|TB. 

The simplest case is the decay rate of the first band Ti, the rate per unit volume 
and per unit time. Since the unit volume contains particles, the one particle rate 
having the dimension of the inverse time is H 



4m 2 
— 9 — 1 

n nig 



This factor 



£i 



32vr \ 



4m 2 



mi 



This exactly coincides with the decay rate of one ^ particle computed in the conven- 
tional way using the Yukawa coupling to (p, | Qym^ ^p 2 . 

The next thing to be checked is the decay rate from the 2nd instability band 
T 2 oc £q- This time one divides the quantity T 2 by n| since two £ particles are 
involved, 



£2 



{9l-9 2 y? 1 



16tt 



m. 



\ 



m c 



m% 



(19) 



In ordinary perturbation theory the amplitude for the 2-body process, — > <pip, 
consists of 3 distinct Feynman amplitudes, the contact term with a single quartic 
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coupling, the t-channel and the u-channel exchange diagrams with two vertices of 
the Yukawa coupling, adding to 

- 2ig\ - i^mK—Lj + — ^) . (20) 
5 t — m A u — wr 

In the zero momentum limit of — > m^, t — m 2 = u — m 2 — > — m|, giving the total 
amplitude 2i(#y — g|). Working out the phase space factor, one finds out that the 
invariant rate v tc \<j given by flux x cross section is identical to T 2 /nJ given above. 

To proceed to the general n— th order case, one notes that the propagator in the 
zero momentum limit is given by 

(21) 



(Ei=i Pi - <?i) 2 - ™ 2 k(n-k)m\' 

where pi is an initial £ momentum and q\ is one of the final tp momenta. Hence the 
n— th order invariant amplitude containing the Yukawa couplings alone is 

m 2(n-i) 11 k(n-k) m n f 2 [(n-l)\] 2 - 1 ) 

It is not difficult to check that this leads to the invariant rate precisely equal to the 
corresponding decay rate of the n— th band, T n /n^. What is left to be shown is 
then the relative weight of the Yukawa and the quartic contribution. Subdiagrams 
of — ► tptp for the whole n£ — > ipip process contribute with a factor, k ^Zk) ^ or 
the Yukawa coupling case, considering the propagator above, and with —ig\ for the 
contact quartic coupling case. The ratio of these two terms is exactly equal to the 
ratio in 



det 



V 



e =-4k(n-k)[ 1 » Y -9t](— T, (23) 
Ak(n-k) -2ig Y ^ ) n k(n-k) y4jV m c ; 



that appears in the decomposition of the submatrix of C n , eq. |15| for the formula T n . 
This proves our assertion that the decay rate of the n— th band T n /n^ is equal to the 
zero momentum limit of the invariant rate computed in the ordinary perturbation 
theory. 

Needless to say, this interpretation is valid only in the small amplitude limit. 
In applications to realistic problems analytic formula in the large amplitude regime 
is indispensable. In cosmo logical application the parameter region along h = 2# 4 , 
or more precisely h — 2# 4 (= 4 fc2 +^ 2 ) <^ Q 4 with the large amplitude of 6*4 3> 1 
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is most important, whose case has been worked out in ref ||. In the present note 
we clarified the physical meaning of the n— th band decay formula when it is taken 
to the £o^O limit. In view of our result here the large amplitude formula in ref 
H may be considered non-perturbative effect which is directly joined to the present 
perturbative result. 
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